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Abstract. In the frame of indirect dark matter searches we investigate the flux 
of high-energy 7-ray photons produced by annihilation of dark matter in caustics 
within our Galaxy under the hypothesis that the bulk of dark matter is composed 
of the lightest supersymmetric particles. Unfortunately, the detection of the caustics 
annihilation signal with currently available instruments is rather challenging. Indeed, 
with realistic assumptions concerning particle physics and cosmology, the 7-ray signal 
from caustics is below the detection threshold of both Cerenkov telescopes and satellite- 
borne experiments. Nevertheless, we find that this signal is more prominent than that 
expected if annihilation only occurs in the smoothed Galactic halo, with the possible 
exception of a ~ 15° circle around the Galactic center if the mass density profile of our 
Galaxy exhibits a sharp cusp there. We show that the angular distribution of this 7- 
ray flux changes significantly if DM annihilation preferentially occurs within virialized 
sub-halos populating our Galaxy rather than in caustics. 



PACS numbers: 95.35. +d,98.35.Gi,98.35.Jk,98.62.Gq,11.30.Pb,98.56.-p,98.70.Rz 
1. Introduction 

Unveiling the nature and distribution of the dark matter (DM) is a fundamental goal 
of both theoretical and experimental physics today. The latest observations PP allow 
~ 23% of the matter in the universe to be of unknown, non-baryonic nature. One of 
the strongest candidates for the DM is a weakly interacting massive particle that was 
non-relativistic at the epoch of decoupling from the rest of the universe, as required by 
the currently popular cold dark matter (CDM) cosmological model. A good DM candi- 
date is represented by the lightest supersymmetric particle (LSP) that in most particle 
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physics scenarios is the neutralino x- This particle is a spin-| Majorana fermion con- 
stituted of a linear combination of the neutral gauge bosons and neutral Higgs doublet 
spartners, \ = aB + bW 3 + cH% + dH%. In the popular SUGRA or SUGRA-like mod- 
els, where gaugino-universality is required, its mass is constrained by LEP limits to be 
greater than ~ 50 GeV j2] while theoretical considerations place an upper limit of a 
few hundred GeV - few TeV [HI H] to the neutralino mass. If i?-parity is conserved, 
neutralinos can change their cosmological abundance only through annihilation. Hence, 
high density regions constitute preferential locations for indirect DM searches. 

Most of the recent works have assumed that the bulk of DM is constituted of neu- 
tralinos and looked for a possible detection of 7-ray annihilation flux with appropriate 
detectors. Phenomenological approaches to DM searches have focused mainly on the 
Galactic Center (GC) that has been observed by the EGRET satellite in the energy 
range 100 MeV - 10 GeV 13 and, more recently, by the ground-based Cerenkov tele- 
scopes WHIPPLE [6 , CANGAROO-II [7] and HESS with an energy threshold of 
hundreds of GeV. These data have shown an enhanced 7-ray flux along the the GC 
direction that has been variously interpreted and, although puzzling, does not yet con- 
stitute a firm evidence of indirect DM detection. New data from VERITAS [5] and 
MAGIC ^0] Cerenkov telescopes, as well as from the soon-to-come satellite experiment 
GLAST [TT], will certainly help clarifying this issue. 

If dark relics are stable and their annihilation cross-section is non-negligible, then 
regions of enhanced density constitute preferential locations for the annihilation of DM 
particles. Our aim is to quantify the DM annihilation rate and to predict the associated 
7-ray flux occurring within those singularities that arise in irrotational flows in a cold, 
collisionless medium during the process of structure formation, usually referred to as 
"caustics" because of their analogy with caustics in geometrical optics ^2] • Although 
the presence of DM caustics is a general prediction of CDM cosmologies, model details 
depend on the characteristics of the DM candidate and on our ability in following the 
build up of caustics in the very nonlinear dynamical regime fBl 113 • In this work 
we will assume that the DM is composed of neutralinos and that the shape and position 
of Galactic caustics can be modeled as in [TH] . 

The rationale behind this work is to predict the annihilation signal from DM caus- 
tics within our Galaxy, characterize its spatial signature and evaluate whether currently 
available or next-generation 7-ray detectors will be able to detect this signal and dis- 
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criminate it from the 7-ray Galactic background. Our results will be compared with 
those obtained from a previous similar analysis performed by [T71 [""J""] that, however, 
considered the case of a single, nearby caustic of non-negligible transverse size. 

The centers of virialized DM halos also constitute possible sites for annihilation of 
DM particles if a singularity in the mass density field called "cusp' develops there as a 
result of the gravitational instability process. Theoretical considerations and numerical 
experiments have neither been able to set strong constraints on the shape of the density 
profile near the halo center [T""J 1201 EU 1221 E3 123 EH] nor to predict the abundance of 
a population of dark galactic sub-halos and their survival time within the Galactic halo 
|2fi | 12 7 | [*"""" ( [*""") ] I3()j . On the other hand, current observations such as the rotation curves 
in the central galactic regions "3~"f] I?""""! 13*3*] . Planetary Nebulae kinematics [31], microlens- 
ing events toward the Galactic bulge [HE] and spectroscopic gravitational lensing jHH] do 
not allow yet to discriminate among the various scenarios, although they seem to favor 
a very shallow, or cored density profile. Nevertheless, many authors have considered 
the possibility of detecting the products of DM annihilation occurring within Galactic, 
extragalactic and sub-galactic DM halos, taking into account both observational and 
theoretical constraints [33 EE1 123 EH H3 El 112] ■ 

As neither theoretical arguments nor numerical experiments have been able yet 
to construct self consistent dynamical models that account for both DM caustics and 
virialized substructures populating a host DM matter halo it is difficult to quantitative 
compare the annihilation flux produced within either structures. Nevertheless, as we 
discuss in this work, qualitative differences between the angular distribution of the sig- 
nals exist that may help in discriminating among these two scenarios. 

The layout of the paper is as follows. In Section 2 we describe the expected 
photon flux from neutralino annihilation and introduce the detectability parameters 
for Cerenkov telescopes and satellite-borne experiments. In Section 3 we study the 
effect of a population of DM caustics on the expected 7-ray flux. Section 4 is devoted to 
investigating the possibility of detecting high energy photons produced in DM caustics. 
In Section 5 we summarize and discuss our main conclusions. 
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Following the prescriptions given in jHZj, we factorize the diffuse photon flux from 
neutralino annihilation as 

^(E 7 , ^, An) = ~ 1[E ^{E 1 ) x $™(^, AO), (1) 

where ip is the angle of view from the GC and AQ is the solid angle of observation for 
the detector. We will not consider here the 7-line emission. Its branching ratio is in 
fact as small as 10 -3 and it produces a very small photon flux which we will neglect. In 
Eq. ^the particle physics is embedded in the term: 

where (cr ann t> ) is the thermally-averaged neutralino self-annihilation cross-section times 
the relative velocity of the two annihilating particles, dn^/dE-y is the differential photon 
spectrum for a given annihilation final state /, whose analytical expression is given in 
|37j . bf is the branching ratio of annihilation into / and m x is the neutralino mass. 
Astrophysics and cosmology affect the last factor of Eq. that represents the line-of- 
sight integral: 

$™(^,AO)=/ dSIf f p 2 x (r(X^'))d\(r^'), (3) 
J An Acs 

where p x (r) is the neutralino density profile, r is the distance from the GC and is equal 
to r = \J X 2 + Rq 2 — 2XRq cosip, X is the distance from the observer and R Q is the 
distance of the Sun from the GC. 

For distant, extragalactic objects, Eq. El reads 

1 /•min[R G ,r max (An)] 

$™(^, AQ) = -j 47rr 2 p 2 x (r)dr, (4) 

where d is the distance of the object from the observer, R G is its virial radius and 
r max (Af2) is the radius seen within the angular acceptance of the detector. 

As mentioned before, the mass distribution near the center of the DM halo is poorly 
constrained. High resolution numerical experiments have shown that a cuspy density 
profile develops near the halo centers with an inner slope that is bracketed between the 
NFW profile (hereafter NFW97) [IH] 

„NFW97 

NFW97 _ A^s /k\ 

f„Ar-NFW97Wl j_ W<r,NFW97^ 2 ^ ' 



( r / r NFW97) (l + r / r NFW97) 

and the steeper Moore et al. profile (M99) [20 : 
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Scale quantities are fixed by the observations (virial mass or peak rotation velocity of the 
halo) and by the estimation of the concentration parameter c = r vir /r s , where we define 
the virial radius r w - 1T as the radius within which the halo average density is 200 times the 
critical density of the universe. The parameters cnfw97 an d CM99 = 0.64 cnfw97 have 
been computed according to [3H] assuming a CDM universe. 

As discussed in Section 1, the existence of a central cusp has not been confirmed by 
recent astronomical observations that, instead, seem to favor the case of a central core 
of constant density. Cored density profile would significantly depress the annihilation 
rate in the central region |44j . However, since using a cored rather than a cuspy density 
profile would not significantly affect our flux predictions for ip > 15°, where the caustics 
signal is expected to be more prominent, in this work we only consider the M99 and 
NFW97 model density profiles. 

Figffl shows the factor $ cosmo (^) of Eq. El for our Galaxy and the three more 
prominent galaxies within a distance of 1.5 Mpc, computed for a solid angle of 10~ 5 sr 
(left panel) and 10 -3 sr (right panel) assuming both a NFW (filled dots) and a Moore 
profile (filled triangles). We would like to remark here that the use of a cored profile 
would substancially reduce the expected fluxes. 

Divergences in the l.o.s. integral are avoided by forcing a constant density core of radius 
10 _8 kpc corresponding to the distance at which the self-annihilation rate equals the 
dynamical time. Several astrophysical effects can increase the size of this core radius, 
leading to a significant decrease in the annihilation rate and expected flux |4TH 13 7j. As 
already pointed out by [40J, these external galaxies shine above the Galactic foreground 
as clearly seen in both panels of Fig^ 

The "supersymmetric factor" d ^ SY in Eq. |2]has been investigated extensively by 
|37j . In our calculations we assume $ SUSY = 10 -32 cm 3 s _1 GeV -2 sr -1 which represents 
a reasonable though somewhat optimistic choice among the values allowed by super- 
symmetric models. 

2.1. Experimental detectability 

One of the aims of this work is to study the experimental detectability of DM 7- 
ray sources with both Cerenkov telescopes and satellite-borne experiments. Following 
[231301, we explore the detectability of the diffuse photon flux from DM annihilation by 
comparing the number n 1 of expected 7 events from neutralino annihilation with the 
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Figure 1. "Cosmological factor" <J> cosmo (^) for the expected 7-ray flux from the most 
prominent external galaxies and for the Milky Way foreground within a solid angle 
Ail = I0~ 5 sr (left panel) or Ail = 10~ 3 sr (right panel), for a M99 (solid line and 
triangles) and a NFW97 (dotted line and circles) profile. A small central core region 
of radius 10 -8 kpc has been assumed. 
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To estimate the experimental sensitivity, a, we choose the subsequent, rather typical 
set of parameters. For the effective observation time, Tg, we consider a 20-day pointing 
for Cerenkov telescope and a 30-day observation for the satellite-borne experiment (for 
a Cerenkov telescope, only periods corresponding to a source zenith angle 9 < 60° are 
considered). The fraction of signal events e^o, is set equal to 0.7 within the optimal 
solid angle AQ = 10~ 5 sr, corresponding to an angular resolution of the instrument of 
0.1°. The effective detection areas for electromagnetic and hadronic induced showers are 
considered to be independent from the value of the particle energy and incident angle 
and are set equal to ^4^ rcnkov = 4 x 10 8 cm 2 and Satellite = 10 4 cm 2 , respectively. The 
identification efficiency e is of crucial importance to reduce the physical background 
level. A Cerenkov telescope has a typical identification efficiency e e . m . ~ 99% for 
electromagnetic induced (primary 7 or electrons) showers and ehad ~ 99% for hadronic 
showers. A satellite experiment has an identification efficiency for charged particles of 
^charged ~ 99.997% and of e neu trai ~ 90% for photons, due to the backslash of high energy 
photons [13]. Our expected signal is compared to the presence of proton jl^l, electron 
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3. The caustics scenario 

The caustics phenomenon has been extensively studied by [5U]. Here we briefly outline 
the main characteristics of the DM caustics that are relevant for our work. 

Collisionless dark matter particles move in the universe on a thin sheet in the 
velocity-position phase space, whose thickness is determined by their small primordial 
velocity dispersion. For a neutralino such a velocity has been calculated, after the 
equivalence time and before the onset of galaxy formation, to be as small as Sv x (t) ~ 
10- n c (GeV/m x ) 1 / 2 (t /t) 2/3 , where t is the age of the universe. Such a sheet winds up 
where an overdensity grows non-linearly (i.e. when a galaxy forms) and the particles on 
the sheet fall onto the corresponding overdensity. The distance from the center of the 
overdensity at which the folding occurs has been estimated to be of the order of 1 Mpc 
for a galaxy like the Milky Way. 

The following period of continuous cross-infalls lasts until inhomogeneities diffuse 
the process and prevent the formation of structures. In this period, a set of caustic 
surfaces which wrap around the galaxy is formed where the phase space folds. These 
surfaces define areas of enhanced density located at the boundary between regions 
characterized by different numbers of flows of particles. Namely, in the internal region of 
the caustic surface there is an odd number n of streams of particles compared to the n — 2 
flows in the external region. The non-zero velocity dispersion prevents the mass density 
to diverge at the caustics surface. More precisely, the density at the caustics is smeared 
off on distances 5x x ~ 3 ■ 10 15 cm (GeV/m x ) 1//2 , corresponding to about 10 -4 pc for a 
TeV-mass neutralino. The final outcome of this process of caustics formation is the onset 
of a series of caustics with different characteristics. A number of outer caustics form, 
that appear as topological spheres surrounding the galaxy, as well as inner caustics, 
which look like tubes lying in the Galactic plane, composed of the particles which carry 
the most angular momentum when they are at their distance of closest approach to the 
Galactic Center. At the intersection between the Galactic plane and a perpendicular 
plane passing through the GC, these caustic tubes have a triangular shaped cross-section 
with three density cusps at the vertices of the triangle (the so-called tricusp), one on the 
Galactic plane and the other two located at the extremes of a segment perpendicular to 
it. 
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3.1. Density distribution near caustics 

The particle density at the caustics can be calculated analytically ^H]- In proximity 
of the cusps of the triangle the density profile d cau scales like r _1 while near both the 
surface of the outer caustics and the side of the triangle it scales like r _1//2 , with r being 
the distance from the caustics. A set of simplifying assumptions is usually made to 
model the density distribution in these regions. 

The claimed observational evidence of a Galactic tube characterized by a triangle with 
side of 100 pc ^H] is at variance with the theoretical prediction of 1 kpc for the side of the 
triangle of the nearest inner caustics, suggesting that neither observations nor theory 
are currently able to set precise constraints on the size of caustic tubes. Therefore, we 
make the assumption that for all practical purposes the caustic triangles are of negligible 
sizes, i.e that the inner caustics reduce to caustic rings, or lines, on the Galactic plane. 
Such an approximation is valid |TH| as long as the distance of the observer from the 
caustic ring is of the same order as the ring radius. This means that, assuming a 
galactocentric distance of 8.5 kpc, caustic rings with radii smaller than about 4-5 kpc 
can be considered as lines. Yet, the approximation still holds for caustic rings which 
are closer to the observer, provided that their transverse dimensions are smaller than 
the distances from the observer which, in this case, are smaller than the ring radius. In 
the model proposed by [THj the nearest caustics is characterized by a radius of 7.8 kpc 
and transverse dimensions of ~ 200 pc. When compared to the galactocentric distance, 
we have that 200 < 700 < 7.800 pc, i.e. both the above conditions are satisfied and the 
approximation is valid. The previous relation is also valid for the remaining caustics 
that have similar transverse dimension but have larger radii. To compute the profile of 
the 7-ray emission from neutralino annihilation in the caustics we furthermore assume 
that the turnaround sphere is initially rigidly rotating and that the velocity of each 
particle in the flow is constant. Under these hypotheses we obtain a simple expression 
for the density distribution near the caustics. 
For the the n-th outer caustics we have EO] : 



where r is the distance from the GC, t> ro t is the rotational velocity of the galaxy, R n is 
the radius of the outer caustics and dM/dQdt is the rate at which mass falls in per unit 
time and unit solid angle, that can be estimated using the self-similar infall model (see 
for detailed references). 




(8) 
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The density profile near the m-th caustic ring is given by [TH I5U] : 




dM 2 1 



(9) 



d ^dt m v v /(r 2 -a^) 2 + 4a^ 2 



where (p, z, 6) are cylindrical coordinates with respect to the GC, a m is the caustics ring 
radius [IS], r = V p 2 + z 2 and v is the velocity amplitude of the DM particles at the 
caustics position. 

3.2. Expected 7 -ray flux from caustics 

Eqs. |H1 and M can be used to compute the expected 7-ray flux from the caustics in our 
Galaxy. The resulting value of the "cosmological factor" $ cosmo (■?/;) computed at the 
Galactic plane is shown in Fig. |21 for a solid angle of 10 -5 sr (left panel) and 10 -3 sr 
(right panel). The expected flux from caustics is obtained by integrating Eq. El over a 
region in which the mass density, given by either Eqs. |H1 or exceeds that of the 
underlying smooth Galactic halo. The mass within caustics, hence the amplitude of the 
annihilation signal, only depends on the virial mass of the host halo (or its rotational 
velocity). Indeed, in the self-similar infall model used to derive Eqs. IHland El the value 
of v ro t determines both the rate of infall mass and caustics' distances from the center of 
the halo. 

The contribution from the outer caustics is negligible on the scale of the plots. A 
most remarkable spatial modulation signature is present which is due to the presence 
of the inner caustics rings (n > 5) with an associated signal well above the Galactic 
foreground level for both the M99 (straight line) and the NFW97 (dashed line) model 
density profiles. 

The contribution of the caustics to the annihilation signal is negligible for 
extragalactic objects. In Fig. |21 we show the case of M31 (filled triangle and starred 
symbol), for which the effect of concentrating part of the halo mass in caustics actually 
decreases the flux. This derives from having placed a sizable fraction of the M31 mass 
in the outer caustics falling outside the viewing angle of the instrument. 

Fig. El shows the ratio (<fi MW + <p cau ) /(f) MW , which defines the enhancement of the 
7-ray signal due to the presence of the caustics rings, for a M99 (straight line) and a 
NFW97 (dotted line) profile. As expected, no enhancement is found in the direction of 
the GC, where the halo emissivity is very large, while DM annihilation in the nearest 
caustics ring results in an enhancement factor of about 600 along the direction ip = 60°. 
The spatial signature of the 7-ray flux expected in presence of Galactic caustics is indeed 
rather prominent. 
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Figure 2. "Cosmological factor" for the expected 7-ray emission from neutralino 
annihilation in the Galactic plane due to the presence of caustics rings for an 
angular resolution of 0.1 degrees corresponding to 10 -5 sr (left panel) and 1 degree 
corresponding to 10~ 3 sr (right panel). The effect of the presence of caustics in the 
M31 galaxy is shown as well (symbols). 




Figure 3. Enhancement of the 7-ray emission from neutralino annihilation in the 
Galactic plane due to the presence of caustics, for a solid angle of 10~ 5 sr. 
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3.3. Ring vs. Tricusp 

As we have anticipated, neither theory nor observations are able to set strong constraints 
on the actual size of the caustics that ranges from ~ 100 pc to ~ 1.0 kpc. It is therefore 
worth investigating how our model predictions change when dropping our hypothesis of 
a ring caustics with negligible transverse siize. This is precisely the case studied by [Tj\ 
that considered the case of the nearest inner caustics (with n=5) located at 7.8 kpc 
from the Galactic Center with transverse dimensions of 1.2 kpc. With this particular 
geometry the Earth, placed at 8.5 kpc from the GC, is located well inside the tricusp. 

The model [T7| also assumes that m x = 100 GeV instead of our fiducial value of 
1 TeV, resulting in a larger smearing distance Sx x and in a smaller cut-off on the caustic 
density profile of ~ 800 GeV cm -3 . Decreasing the density cutoff reduces the expected 
annihilation flux. The reduction factor weakly depends on the geometry of the problem. 
In the case of a sigle caustic ring with negligible size, he value of $^ s ^ decreases by a 
factor of ~ 2, where $^ s ^ is the maximum of $ cosmo (■?/>, Afi = 10~ 5 sr _1 ). In the case 
of a tricusp geometry considered by [T7] the effect is very similar and decreases 
by a factor of ~ 3.3. 

Varying the geometry while keeping m x constant has a much more dramatic im- 
pact. We have found that with m x = 100 GeV the value of $^ s ^ for a tricusp is ~ 100 
smaller than that expected from a caustic ring with negligible transverse size. This 
effect is hardly surprising as in the tricusp scenario, and unlike in the ring case, only a 
small portion of the caustics is seen within the angle of view AQ. 

Overall, the value of ^^amo for the caustics n = 5 computed as in Section 3.2 is 
expected to be ~ 100 times larger than that of [Hj. Note, however, that discrepancy in 
the flux predictions of the two models are less severe. This derives from the fact that 
our value of $ SUSY = 10~ 32 cm 3 s" 1 GeV -2 sr _1 is ~ 1000 smaller than that assumed by 
[TTj which, however, requires very optimistic hypotheses about the underlying particle 
physics. 

4. Possibility of 7-ray flux detection 

To investigate the possibility of detecting the flux of 7-ray photons produced by 
neutralinos annihilating in DM caustics, we compute the sensitivity of both a Cerenkov 
ground-based telescope and a satellite-borne experiment to the annihilation fluxes shown 
in Fig. 121 The parameters used to compute the experimental sensitivity (Eq. EJ) are 
specified in Section |2~T1 and represent realistic observations that could be carried out by 
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next-generation detectors. 

Fig. 0] shows the experimental sensitivity, a, a function of the angle from the GC 
n the case of a Cerenkov telescope at the VERITAS site (left panel) and a GLAST-like 
satellite (right panel). The choice of the Cerenkov telescope site determines the observ- 
able region of the sky and thus may constitutes a source For instance the GC region, 
which is maybe the brightest spot in the 7-ray sky, cannot be observed at the VERITAS 
site considered in this work. It is worth stressing that, because of its narrow field of 
view, each point in the left panels of Fig. H] represents a single, 20-day long observation 
with a Cerenkov telescope. In both panels the significance of the detection is is pro- 
portional to the 7-ray flux and thus exhibits the same dependence from found in Fig.EJ 

Clearly, the significance detection level is very low, even in correspondence of the 
peaks of the experimental sensitivity. This means that 7-ray photons produced by DM 
annihilation in the caustics cannot be detected by currently planned experiments unless 
new theoretical input from particle physics could boost up the supersymmetric factor 
by a factor of 30, which seems rather implausible. 

These pessimistic conclusions are not alleviated when taking into account for the 
large theoretical uncertainties in the caustics model. Indeed, as we have shown in Section 
3.3, increasing the transverse dimension of the nearest caustics while keeping constant 
its distance from the GC and using the same value of m x , has the effect of reducing 
the expected annihilation signal by a factor of ~ 100 making, it more challenging its 
experimental detection. Moreover, a further decrease of the annihilation signal is to 
be expected if dynamical disturbances within the Galactic halo are large enough to 
significantly shorten their survival time (|27j). 

5. Conclusions 

In the popular CDM cosmological model, high density sub-structures are expected to 
form via gravitational instability from small perturbations in the primordial density field 
of cold, collisionless particles. In this work we have investigated the annihilation of DM 
particles within caustics, that represent regions of very high density, and have estimated 
the possibility of detecting the 7-ray photons produced by annihilation events within our 
Galaxy. Our main assumptions concern the geometry of the inner caustics, that we have 
considered small enough to be described as a series of line rings lying in the Galactic 
plane at various distances from the GC, as well as the underlying particle physics, namely 
the value of the neutralino mass m x = 1 TeV and of $ SUSY = 10 -32 cm 3 s _1 GeV~ 2 sr _1 . 
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Figure 4. Left panel: expected sensitivity of a Cerenkov telescope to 7-ray emission 
from neutralino annihilation in the Galactic plane due to the presence of caustic rings, 
at an energy threshold of 100 GeV for a solid angle of 10~ 5 sr. Right panel: the same 
as in the left panel calculated for a satellite-borne experiment and an energy threshold 
of 10 GeV. 



Our main conclusion is that the expected annihilation signal from inner Galaxy 
caustics cannot be revealed with presently available or currently planned Cerenkov 
telescopes and satellite-borne detectors. This result takes into account the present 
uncertainties in modeling the caustics geometry and their spatial distribution, both of 
which systematically contribute in suppressing the expected annihilation signal. As 
varying m x does not change significantly the value of $ cosmo ; the only possibility of 
observing an annihilation signal from DM caustics would be that of assuming a much 
larger value for $ SUSY . For example, ^] found that a significant detection is 
indeed possible in one year of observation with GLAST if a very optimistic value 
of $ SUSY = 10~ 29 cm 3 s~ x GeV -2 sr _1 is assumed that, however, we regard as rather 
implausible. 

Although the possibility of its detection is small, some additional considerations 
about the annihilation signal are worth being made. First of all, the annihilation flux 
from Galaxy caustics dominates over the smooth Galactic background for ip > 15° and 
possibly also within this region if a constant density core, rather then a sharp cusp, 
is present in the GC. Moreover, as we have checked, the annihilation signal from DM 
caustics does not significantly contribute to the total annihilation signal expected from 
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extragalactic sources such as M31. 

Finally, it is interesting to compare the annihilation flux in the caustics scenario 
with that expected from a population of virialized sub-galactic halos that has already 
been studied by 123 I2H EIE EH] • Indeed, the annihilation signal produced in Galaxy 
caustics is expected to have an angular distribution across the sky that is very different 
from that produced in the smooth Galactic halo (Fig. EJ) and from that emitted by a 
population of virialized sub-galactic halos. In Fig. E] we show the enhancement factor 
within a solid angle of 10 -5 sr due to the presence of sub-Galactic halos of masses 
> 10 6 M Q containing 20% of the original halo mass and having the same NFW97 density 
profile as the Galactic halo. The enhancement factor obviously depends on the modeling 
of the sub-Galactic halos density profile, mass function and spatial distribution (e.g. the 
three histograms of Fig. Elshow the effect of using the different models of sub-halos tidal 
disruptions proposed by |4*0]). Moreover, we would like to point out again that the use of 
a cored subhalos' density profile instead of the proposed NFW97 would further decrease 
the expected fluxes. Yet, as we have verified, accounting for model uncertainties does 
not significantly affect the angular distribution of the annihilation signal, in the sense 
that the enhancement profile has a peak that is broader than that predicted by the 
caustics scenario with a maximum at larger angles from the GC (ip ~ 120° instead of 
~ 60°). The enhancement level remains low too, which means that 7-ray emission from 
sub-Galactic halos will hardly be detected by future gamma-ray detectors, especially 
because observations by satellite-borne experiments like GLAST will be hampered by 
the very long exposure time required by the serendipitous nature of these 7-ray sources. 

Based on these considerations, the GC remains the best place to look at for 
detecting 7-ray signal, as it is shown in Fig. IH provided that a sharp DM cusp is 
present there. However, ground-based observatories located in the northern hemisphere 
could only observe the GC at a large zenith angle, hence introducing large observational 
errors. For Cerenkov telescopes located in the southern hemisphere, such as HESS 
and CANGAROO, as well as for satellite-borne detectors, the GC represents by far 
the brightest source of the 7-ray sky and will therefore give the best chance to detect 
a possible signature of neutralino annihilation. Much more sensitive detectors will be 
needed to investigate the nature of the DM sources responsible for such signal. 
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